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ABSTRACT 

We have obtained HST FOC f/48 long-slit spectroscopy of the inner 4 arcseconds of 
the Narrow Line Region of NGC 1068 between 3500-5400A with a spectral resolution 
of 1.78A/pixel. At a spatial scale of 0'.'0287 per pixel these data provide an order of 
magnitude improvement in resolution over previous ground based spectra and allow 
us to trace the interaction between the radio jet and the gas in the NLR. Our results 
show that, within ±0'/5 of the radio-jet the emission lines are kinematically disturbed 
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and split into two components whose velocity separation is 1500 km s _1 . The filaments 
associated with the radio lobe also show a redshifted kinematic disturbance of the order 
of 300 km s _1 which probably is a consequence of the expansion of the radio plasma. 

Furthermore, the material enveloping the radio-jet is in a much higher ionization 
state than that of the surrounding NLR gas. The highest excitation is coincident with 
the jet axis where emission in the coronal line of [FeVII] A3769A is detected and the 
HeIIA4686 A is strong but where [Oil] A3727 A is depressed. This large localized 
increase in ionization on the jet axis is accompanied by the presence of an excess con- 
tinuum. Because the electron density is substantially larger in the jet compared to the 
surrounding NLR, these results can only be explained if there is a more intense ionizing 
continuum associated with the jet. This can be accomplished in a variety of ways which 
include an intrinsically anisotropic nuclear radiation field, a reduced gas covering factor 
or the presence of a local ionization source. 

The morphology, kinematics and, possibly, the ionization structure of the NLR in 
the vicinity of the jet of NGC 1068 are a direct consequence of the interaction with the 
radio outflow. 

Subject headings: Galaxies - individual (NGC1068); Galaxies - Seyfert; Galaxies - active 
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1. Introduction 

Extensive HST emission line imagery of the 
Narrow Line Regions (NLR) of Seyfert galaxies 
with linear radio sources has shown that the mor- 
phology of the NLR is directly related to that 
of the radio emission. Seyfert galaxies with a 
lobe-like radio morphology (e.g. Mrk 573, Mrk 
78, NGC 3393, IRAS 0421+045 and IRAS 1105- 
035) have bow-shock shaped emission line regions 
(|Bower et al. 1994j , |Axon et al. 1996| , |Capetti et 
|al. 1996| ), while those with a jet-like radio struc- 
ture (e.g. Mrk 3, Mrk 348, Mrk 6, Mrkl066) 



have jet-like emission line struct ures (|Bowcr et 
al. 1995 , |Capetti et al. 1995a , Capetti et al 



1995b| ). These results have provided compelling 



evidence for strong dynamical interactions be- 
tween the emission line gas and radio-emitting 
ejecta on scales of < 1 kpc. 

The physical foundation of this picture was 
first developed by Taylor, Dyson and Axon (1992) 
who constructed fast bow shock models for the 
interaction between radio ejecta and the ISM 
which included the effects of photo-ionization 



from the nucleus. Dopita and Sutherland ( Do- 



emitted by undisturbed gas in the disk. At larger 
radii this triple structure disappears. 

The complexity of the velocity field observed 
in the inner region is almost certainly due to the 
spatial confusion caused by projection and seeing 
effects which mixes kinematic components from 
different regions along the line of sight. In par- 
ticular the velocity structure created by the ex- 
pansion of the hot gas associated with the lobe is 
projected onto the region interacting with the jet 
which can be spatially resolved only with HST. In 
order to obtain a more detailed picture of the re- 
lationship between the velocity field and the ion- 
ization conditions in the NLR gas and the inter- 
action of the gas with the radio jet we have there- 
fore obtained new long-slit spectroscopic obser- 
vations with the Faint Object Camera f/48 spec- 
trograph on board the HST. Our results allow us 
to resolve the spatial confusion and provide evi- 
dence that shocks created by the jet play a key 
role in forming and, possibly, ionizing the NLR of 
NGC 1068. Throughout this paper we will adopt 
a distance to NGC 1068 of 14.4 Mpc ( [Tully 1988]) 
where 1" corresponds to 72 pc. 



hita fc Sutherland 19964 Popita fc Sutherland! 2 . Observations and Data Reduction 



1996bj ) have emphasized that the hot shocked gas 



could itself make an important local contribution 
to the ionizing continuum (becoming what they 
termed an auto-ionizing shock). 

NGC 1068 is one of the closest Seyfert galax- 
ies and it harbours a bright radio source, with 
a prominent radio-jet ( Muxlow et al. 1996j ) ter- 



minating in an extended radio- lobe (Wilson & 



Ulvestad 1983[ ). Ground-based kinematic stud- 
ies of the NLR have been carried out previously 



by a number of authors ( Alloin et al. 1983 , Bald- 



|w 
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in et al. 1987, Meaburn and Pedlar 1986, Ce- 



il, Bland | Tully 199fl , |Unger et al. 1992- , |a7 
|ribas et al. 19961) in [OIII] A5007. In particular 
two broad components straddling a narrow cen- 
tral component in the NE radio lobe were iden- 
tified. The narrow component appears to follow 
the expected rotation curve of the galaxy and is 



The Narrow Line Region of NGC 1068 was 
observed using the FOC f/48 long-slit spectro- 
graph on October 18 th , 1996 at resolutions of 
1.78A and 070287 per pixel along the dispersion 
and slit directions, respectively. The F305LP 
filter was used to isolate the first order spec- 
trum which covers the 3650-5470 A region and 
therefore includes the [OII]A3727, H/3A4861 and 
[OIII]AA4959, 5007 A emission lines. The slit, 
/ /063xl3'. / 5 in size, was placed at a position an- 
gle of 83° and spectra with exposure times of 627 
seconds were taken in the 1024x512 non-zoomed 
mode at 6 locations separated by 1". Unfortu- 
nately, due to a guide star re-acquisition failure, 
only 3 of the slit locations yielded usable spec- 
tra which we identify as POS1, POS2 and POS3 
respectively. 
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Fig. 1. — The slit positions used in our study overlayed on a gray scale map of the [OIII] image of 
Macchetto et al. 1994. The origin of the spatial axis on the borders of the image is centered on the 
position of the hidden nucleus (Capetti et al. 1995c). The spatial coordinates used along the slit are 
marked as crosses at POS1 and POS2. The zero point of this axis lies on cloud G in POS3. The filaments 
associated with the radio lobe intersect the slit at POS2 at location -1.5 arcseconds and POS1 at -1 
arcseconds. The filled squares at POS3 indicate the limits of the contiguous regions where we extracted 
the spectra shown in Fig. 3. 
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The data reduction follows the procedure de- 
scribed in detail by Macchetto et al. 1997 , Flux 



calibration was performed using the observa- 
tions of the spectrophotometric standard star 
LDS749b. 

To accurately determine the location of the 
slits we compared the surface brightness profile 
derived from the FOC, f/96, F501N HST image 
of Macchetto et al. 1994] with that measured from 
the spectra at the three slit positions. The best 
match is displayed in Fig. [j] and is accurate to 
within half a slit width (~ 0"03). 

The slit locations POS1 and POS2 cut across 
the emission line filament system associated with 
the Northern radio lobe ( papetti, Axon, &z Mao 



che jtto 1997| ) . While the slit at POS3 (see Fig. 
H) traverses cloud G ( Evans et al. 1991~| ) and the 
parallel feature to the West which form a funnel- 
like structure where the jet meets the radio lobe 



(see Gallimore et al. 1996). 



3. Results 

The velocities derived from the [OIII] and 
[O II] lines by fitting gaussian line profiles are 
plotted in Fig. |2[ Around the jet axis, at POS3 
(Fig. ||b), the lines are split into two velocity sys- 
tems separated by ~ 1500 km s _1 . Outside this 
region the gas motions are quiescent but with a 
gradient of +500 km s^ 1 from West to East. No- 
tice that the velocity splitting is not symmetric 
about the systemic velocity, but is much larger on 
the blue-shifted side. The blueshifted component 
is also substantially brighter and has a larger 
velocity dispersion than the redshifted compo- 
nent. At POS2 (Fig. |a) the velocity field is 
relatively flat, except at the filaments at slit co- 
ordinates ~ -1" - -2", where there is a large red- 
shifted velocity perturbation of amplitude ~ 300 
km s -1 , at a similar radial velocity to that of the 
redshifted component of the split line region in 
POS3. POS1 (Fig. ||a) shows a similar pattern 
in velocity, but with a reduced amplitude in the 
perturbed region. 



Turning to the emission line strengths, in Fig. 
H we show plots of representative spectra from 
the five regions identified in Fig. H. On the jet 
axis (labeled as "Jet") the gas is highly ionized; 
[Fe VII]A3759A and [Ar IV] are clearly visible and 
the HeIIA4686 is as high as 0.7 the strength of 
H/3, while the [Oil] emission appears to be de- 
pressed. This increase in excitation is accompa- 
nied by the presence of an excess local contin- 
uum, as shown in Fig. ||. 

Outside the jet the excitation conditions are 
much lower, with [Oil] being relatively stronger, 
particularly to the E as shown in Fig. ^. The 
variation of the [OIII]/[OII] ratio along the slit, 
plotted in Fig. [|, shows that there is a fac- 
tor ~5 difference between the ratio on the jet 
axis and outside. The electron temperature 
measured from the [O III]AA4363/5007 ratio give 
T e = (1.15 ± 0.15) xlO 4 K. A lower limit for 
the electron density can be obtained from the 
[Ar IV]AA4710/4740 ratio and yields N e > 10 4 - 5 cm~ 3 . 
At similar distances from the nucleus, but off 
the jet, ground based density measurements from 
[S II]A6716,6731 ratio, obtained from our unpub- 
lished spectroscopy, indicate electron densities of 
a few 100 cm -3 . Thus the increase in excitation 
on the jet occurs despite an increase in density. 

4. Discussion 

The morphological connection between the op- 
tical and radio emission in the region of the jet 
has been discussed in detail by Gallimore et al. 

(1996) and by Capetti, Macchetto and Lattanzi 

(1997) . They showed that there is a clear anti- 
correlation between radio and optical emission, 
with the brightest line emitting knots surround- 
ing the radio-jet. In the jet-cloud interaction 
picture, the radio jet clears a channel through 
the cool ISM, leaving hot ionized gas in its wake 
and the line emission is highly enhanced along 
the edges of the radio-jet where the gas is com- 
pressed. 

The overall velocity structure and ionization 
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Fig. 2. — Velocity fields measured at the three slit positions. Lower panel: Around the jet axis, at POS3, 
the gas is strongly kinematically perturbed and shows split lines with components separated by ~ 1500 
km s _1 . The dashed line marks the systemic velocity of NGC 1068 (Baan & Haschick 1983). Upper 
panel: at POS2 and POS1 there are also large velocity perturbation at the location of the emission line 
filaments. 
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Fig. 3. — Spectra extracted at the positions indicated in Fig. 1. The spectrum extracted in a / . / 14xO'. / 063 
aperture around cloud G, labeled "Jet", shows the high excitation gas on the jet with [FeVII], [ArlV] 
and He II emission lines. 
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Fig. 4. — Grey scale image of the spectrum at P0S3 showing the excess continuum associated with the 
radio jet. 
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Fig. 5. — A plot of the [O III] /[Oil] ratio as a function of position along the slit at POS 3. The ratio 
reaches a maximum on the jet axis where it is ~ 5 times larger than that in the surrounding regions. 
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conditions of the gas associated with the jet of 
NGC1068 can be understood in terms of this pic- 
ture of an expanding and cooling cocoon around 
the jet which creates very fast shocks with at 
least V — 700 km s _1 (even when mass-loading 
of the expanding cocoon is taken into account). 

The viewing angle is such that to the NE the 
lobe is projected on to the disk of NGC1068, so 
that the jet is seen emerging towards us out of 
the disk. The situation is very similar to that 
seen previously on the jet of 3C120 ( [Axon et 



al. 1989 ), even down to the clumpy structure of 
the gas, which almost certainly arises because of 
the onset of instabilities caused by the mixing of 
hot and cold material. This morphology is again 
readily understood if it is due to the interaction 
between the jet and the surrounding medium. A 
viable explanation for the observed asymmetry of 
the velocity splitting is that it is due to density 
stratification of the disk, that is the expanding 
hot bubble is able to expand more rapidly along 



is approximately proportional to U; Penston et 
al 1989). The enhancements in both density and 
ionization parameter imply that the local ioniz- 
ing radiation is a factor ~ 500 higher at the lo- 
cation of the jet in POS 3, than it is elsewhere. 
The exceptionally strong [Fe VII] and Hell lines 
suggest that matter-bounded clouds might con- 
tribute significantly to the line emission near the 
jet. In this case, we will be overestimating the 
jump in U since low ionization lines such as [Oil] 
are strongly suppressed in such clouds. Neverthe- 
less, this is unlikely to explain the whole of the 
inferred increase in the ionizing flux. Either a lo- 
cal ionization source is present which dominates 
the AGN radiation field near the jet axis, or the 
AGN radiation field is itself highly anisotropic. 

The most important local ionization source 
is likely to be ionizing radiation (free-free and 
high ionization line emission) from gas heated 
and ionized in fast shocks driven by the radio 
jet ( Binette, Dopita fe Tuohy, 1985| ; Sutherland, 



the density gradient away from the disk plane, Bicknell fe Dopita 1993| ; Dopita & Sutherland 



than towards it. Further weight to this argu- 
ment is provided by the systematic blue-shift of 
coronal lines with respect to the lower excitation 



lines ( Marconi et al. 1996 ). This implies that the 
coronal lines comes from a different kinematic 
component which we would ascribe to the hot 
cooling gas in the cocoon. 

The spectral diagnostics obtained from our 
long-slit spectra indicate that the kinematically 
disturbed line-emitting gas associated with the 
radio jet has both a higher density and a higher 
excitation level than the surrounding NLR. A 
similar result, based on narrow band imaging, 



has been reported recently by Capetti, Axon, fe 



M qcchetto 1997| . The density enhancement is a 
factor > 100 if the ground-based measurements 
of the [SII] AA6716,6731 ratio can be taken as typ- 
ical of the regions extending away from the jet 
axis. On the other hand, the variations in the 
[OIII]5007/[OII]3727 ratio along the POS 3 slit 
imply a factor 5 enhancement in the ionization 
parameter (U) across the jet axis (since this ratio 



1996a, b). The strong [Fe VII] emission and the 
excess continuum observed on the jet axis may be 
a clue that the coronal gas expected in this pic- 
ture is indeed present. Alternatively, the on-axis 
increase in the ionizing flux could be attributed 
to intrinsic anisotropy of the nuclear continuum. 
In this case, the size of the enhancement in flux, 
and the fact that it is highly localized, require 
a continuum beaming factor comparable with 
Doppler boosting in a moderately relativistic jet 
(gamma a few). This idea can be ruled out 
because the region we have studied occurs after 
a significant bend of the jet and therefore one 
would expect to see higher excitation along the 
orginal jet direction, which is not observed (cf. 
Capetti, Axon, fe Macchetto 1997| ). 

Perhaps a more likely explanation is that the 
continuum anisotropy is largely produced by az- 
imuthal variations in the optical depth to ioniz- 
ing photons. In this picture, the radio jet sweeps 
out a channel as it passes through the NLR with 
the result that ionizing radiation from the cen- 
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tral source suffers much less attenuation along its 
axis than in other directions. 

Our results on POS1 and POS2 can also be 
interpreted in terms of the interaction between 
the radio plasma and the ambient gas, but in this 
case it is the expanding lobe material which is 
driving the motion. From our spectroscopic data 
we can say little on the ionization conditions of 
these filaments because only the [O III] lines are 
detected. 

5. Conclusions 

We have presented new HST FOC f/48 long- 
slit spectra of the central 4 arcseconds of the 
NLR of NGC 1068. At a spatial scale of / .'0287 
per pixel these data provide an order of mag- 
nitude improvement in resolution over previous 
ground based spectra and allow us to trace the 
interaction between the radio jet and the gas in 
the NLR. Our results show that, within ±0 / /5 of 
the radio-jet the emission lines are split into two 
components whose velocity separation is ~ 1500 
km s _1 and are clearly the result of the interac- 
tion between the radio jet and the ambient NLR 
material. The filaments associated with the ra- 
dio lobe also show a redshifted kinematic distur- 
bance of the order of 300 km s _1 which probably 
is a consequence of the expansion of the radio 
plasma. 

Our results show that the highest excitation 
gas is physically related to the location of the 
jet and is accompanied by an excess continuum 
providing additional evidence for the interaction 
model. 

The morphology, kinematics and, possibly, the 
ionization structure of the NLR in the vicinity of 
the jet of NGC 1068 are a direct consequence of 
the interaction with the radio outflow. 
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